JLAB-THY-00-09 



Realistic Calculation of the hep Astrophysical Factor 



L.E. Marcucci( a \ R. Schiavilla (a ' b) , M. Viviani^, A. Kievsky^, S. Rosati (c ' d) 
( a ) Department of Physics, Old Dominion University, Norfolk, VA 23529 
(b) Jefferson Lab, Newport News, VA 23606 
&INFN, Sezione di Pisa, 1-56100 Pisa, Italy 
^Department of Physics, University of Pisa, 1-56100 Pisa, Italy 

(February 8, 2008) 



Abstract 

The astrophysical factor for the proton weak capture on 3 He is calculated 
with correlated-hyperspherical-harmonics bound and continuum wave func- 
tions corresponding to a realistic Hamiltonian consisting of the Argonne v is 
two-nucleon and Urbana-IX three-nucleon interactions. The nuclear weak 
charge and current operators have vector and axial-vector components, that 
include one- and many-body terms. All possible multipole transitions con- 
necting any of the p 3 He S- and P-wave channels to the 4 He bound state are 
considered. The 5"-factor at a p 3 He center-of-mass energy of 10 keV, close to 
the Gamow-peak energy, is predicted to be 10.1 x 10 -20 keV b, a factor of 
five larger than the standard-solar-model value. The P-wave transitions are 
found to be important, contributing about 40 % of the calculated S'-factor. 
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Recently, there has been a revival of interest in the reaction 3 He(p,e + z/ e ) 4 He Jl[]. This 
interest has been spurred by the Super-Kamiokande collaboration measurements of the en- 
ergy spectrum of electrons recoiling from scattering with solar neutrinos 0. At energies 
larger than 14 MeV more recoil electrons have been observed than expected on the basis of 
standard-solar-model (SSM) predictions ||. The hep process, as the proton weak capture 
on 3 He is known, is the only source of solar neutrinos with energies larger than 15 MeV- 
their end-point energy is about 19 MeV. This fact has naturally led to questions about the 
reliability of the currently accepted SSM value for the astrophysical factor at zero energy, 
2.3 x 10~ 20 keV b ||. In particular, Bahcall and Krastev have shown that a large en- 
hancement, by a factor in the range 20-30, of the SSM S'-factor value given above would 
essentially fit the observed excess g] of recoiling electrons. 

The theoretical description of the hep process, as well as that of the neutron and pro- 
ton radiative captures on deuteron and 3 He, constitute a challenging problem from the 
standpoint of nuclear few-body theory. Its difficulty can be appreciated by comparing the 
measured values for the cross section of thermal neutron radiative capture on 1 H, 2 H, and 
3 He. Their respective cross sections are: 334.2 ± 0.5 mb ||, 0.508 ± 0.015 mb |J, and 
0.055 ± 0.003 mb |7|]. Thus, in going from A=2 to 4 the cross section has dropped by al- 
most four orders of magnitude. These processes are induced by magnetic-dipole transitions 
between the initial two-cluster state in relative S-wave and the final bound state. In fact, 
the inhibition of the A=3 and 4 captures has been understood for a long time |§: the 3 H 
and 4 He states are approximate eigenstates of the magnetic dipole operator //, and conse- 
quently matrix elements of fi z between nd (n 3 He) and 3 H ( 4 He) vanish (approximately) due 
to orthogonality. This orthogonality argument fails in the case of the deuteron, since then 
jj, z can connect the large S-wave component of the deuteron to an isospin T=l 1 S np state. 

This quasi-orthogonality, while again invalid in the case of the proton weak capture 
on protons |§, is also responsible for inhibiting the hep process. Both these reactions are 
induced by the Gamow- Teller operator, which differs from the (leading) isovector spin part 
of the magnetic dipole operator essentially by an isospin rotation. As a result, the hep weak 
capture and nd, pd and n 3 He radiative captures are extremely sensitive to: i) the D-state 
admixtures generated by tensor interactions, and ii) many-body terms in the electroweak 
current operator. For example, many-body current contributions provide, respectively, 50 % 
and over 90 % of the calculated pd [|10j and n 3 He fi],|IHl cross sections at very low energies. 

In this respect, the hep weak capture is a particularly delicate reaction, for two additional 
reasons: firstly and most importantly, the one- and many-body current contributions are 
comparable in magnitude, but of opposite sign PUl2fl; secondly, many-body axial currents, 
specifically those arising from excitation of A isobars which give the dominant contribution, 
are model dependent fi~2"f . This destructive interference between one- and many-body cur- 



rents also occurs in the n 3 He ( "hen n ) radiative capture flJTTH , with the difference that there 
the leading components of the many-body currents are model independent, and give a much 
larger contribution than that associated with the one-body current. 

The cancellation in the hep process between the one- and two-body matrix elements has 
the effect of enhancing the importance of P-wave capture channels. Indeed, one of the results 
of the present work is that these channels give about 40 % of the S-f actor calculated value. 
That the hep process could proceed as easily through P- as S-wave capture was not not 
sufficiently appreciated JT3] in all earlier studies of this reaction we are aware of, with the 
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exception of Ref . [|IJ , in which Horowitz suggested, on the basis of a very simple one-body 
reaction model, that the 3 Po channel may be important. 

Most of the earlier studies 0,0,|T^] had attempted to relate the matrix element of the 



axial current occurring in the hep capture to that of the electromagnetic current in the 
hen capture, exploiting (approximate) isospin symmetry. This approach led, however, to 
S'-factor values ranging from 3.7 to 57, in units of 10~ 20 keV b. In an attempt to reduce the 
uncertainties in the predicted values for both the radiative and weak capture rates, ab initio 
microscopic calculations of these reactions were performed in the early nineties P, |TT| , |T2 



using variational wave functions corresponding to a realistic Hamiltonian, and a nuclear 
electroweak current consisting of one- and many-body components. These studies showed 
that inferring the hep S-factor from the measured hen cross section can be misleading, 
because of different initial-state interactions in the n 3 He and p 3 He channels, and because of 
the large contributions associated with the two-body components of the electroweak current 
operator, and their destructive interference with the one-body current contributions. 

The significant progress made in the last few years in the modeling of two- and three- 
nucleon interactions and the nuclear weak current, and the description of the bound and 
continuum four-nucleon wave functions, have prompted us to re-examine the hep reaction. 
In the present work we briefly summarize the salient points in the calculation, and report 
our results for the S'-factor in the energy range 0-10 keV. An exhaustive account of this 
study [[15|], however, will be published elsewhere. 



The cross section for the 3 He(p,e + z/ e ) 4 He reaction at a cm. energy E is written as 
"(E)-- 




x i£S k/i ^ |!>|2 IfIf' (1) 

where Am = m + m 3 — m 4 = 19.8 MeV (m, 777,3, and m 4 are the proton, 3 He, and 4 He rest 
masses, respectively), t> re i is the p 3 He relative velocity, and the transition amplitude is given 
by 



G 



(f\H w \i) = r|r(-q; 4 He|£(q)| P ;p 3 He) . (2) 

Here Gy is the Fermi constant, q = Pe + p^, |p;p 3 He) and | — q; 4 He) represent, respectively, 
the p 3 He scattering state with relative momentum p and 4 He bound state recoiling with 
momentum — q, l a is the leptonic weak current, l a = u u ^ a {l — 7s)f e (the lepton spinors are 
normalized as v\v e = u^ v u v = 1), and j°"(q) is the nuclear weak current, j CT (q) = (p(q), j(q))- 
The dependence of the amplitude upon the spin projections of the leptons, proton and 3 He 
has been omitted for ease of presentation. Since the energies of interest are of the order of 
10 keV or less-the Gamow-peak energy is 10.7 keV for the hep reaction-it is convenient to 
expand the p 3 He scattering state into partial waves, and perform a multipole decomposition 
of the nuclear weak charge, p(q), and current, j(q), operators. Standard manipulations lead 
to[0 

\ E E I (f I H w I *)? = Gy L(j T N°\ ( 3 ) 

S e S v S1S3 
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where the lepton tensor L ar is written in terms of electron and neutrino four-velocities as 

L ar 



v ff v T + V V 

e v 1 v e 



g aT v e ■ v u + i e CTar/3 v e Q v !/j/ g, while the nuclear tensor is defined as 



with 



(4) 



S1S3 



W°=°' 3 (q; Sl s 3 ) = *o L5J (q; SiS S )Tf SJ (q) 



LSJ 



■-4 E *-f J (q; *i*s) 



LSJ 
rLSJ 



\Mf J (q) + Ef J (q) 



(5) 



(6) 



where A = ±1 denote spherical components. The functions A A=0i±1 depend upon the direc- 
tion q, and the proton and 3 He spin projections s\ and s 3 [|L5|] (note that the quantization 
axis for the hadronic states is taken along p, the direction of the p 3 He relative momentum), 
while T^ SJ =C L j SJ or L L / J for a=0 or 3. The quantities Cf J , Lf J ', E L j SJ and Mj SJ are 
the reduced matrix elements (RMEs) of the Coulomb (C«J, longitudinal (£«J, transverse 
electric (Eu z ), and transverse magnetic {M^ z ) multipole operators between the initial p 3 He 
state with orbital angular momentum L, channel spin S (S=0,1), and total angular momen- 
tum J, and final 4 He state. The present study includes S- and P-wave capture channels, i.e. 
the 1 So, 3 Si, 3 Po, 1 Pi, 3 Pi, and 3 P2 states, and retains all contributing multipoles connecting 
these states to the J 7r =0 + ground state of 4 He. 

The bound- and scattering-state wave functions are obtained variationally with the 
correlated-hyperspherical-harmonics (CHH) method, developed in Refs. P^JI7| . The nuclear 
Hamiltonian consists of the Argonne v\& two-nucleon 
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[p~8|1 and Urbana-IX three-nucleon 
interactions. This realistic Hamiltonian, denoted as AV18/UIX, reproduces the experimen- 
tal binding energies and charge radii of the trinucleons and 4 He in "exact" Green's function 



Monte Carlo (GFMC) calculations [2D]. The binding energy of 4 He calculated with the CHH 
method [0,0 is within 1 % of that obtained with the GFMC method. The accuracy of 
the CHH method to calculate scattering states has been successfully verified in the case of 
three-nucleon systems, by comparing results for a variety of Nd scattering observables ob- 
tained by a number of groups using different techniques [21]. Studies along similar lines [22| 



to assess the accuracy of the CHH solutions for the four-nucleon continuum have already 
begun. 



The CHH predictions [|I7| for the n H total elastic cross section and coherent scattering 
length have been found to be in excellent agreement with the corresponding experimental 
values. The n 3 H cross section is known over a rather wide energy range, and its extrapolation 
to zero energy is not problematic p3 |. The situation is different for the p 3 He channel, for 
which the singlet and triplet scattering lengths a s and at have been determined from effective 
range extrapolations of data taken above 1 MeV, and are therefore somewhat uncertain, 
a s = (10.8 ± 2.6) fm |4| and a t = (8.1 ± 0.5) fm [|| or (10.2 ± 1.5) fm 0. Nevertheless, 
the CHH results are close to the experimental values above: the AV18/UIX Hamiltonian 
predicts |17| a s = 11.5 fm and a t = 9.13 fm. At low energies (below 4 MeV) p 3 He elastic 



scattering proceeds mostly through S- and P-wave channels, and the CHH predictions, based 
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on the AV18/UIX model, for the differential cross section [^5| are in good agreement with 
the experimental data. 

The nuclear weak current has vector and axial-vector parts, with corresponding one- 
and many-body components. The one-body components have the standard expressions 
obtained from a non-relativistic reduction of the covariant single-nucleon vector and axial- 
vector currents, including terms proportional to 1/m 2 . The two-body weak vector currents 
are constructed from the isovector two-body electromagnetic currents in accordance with the 
conserved- vector-current hypothesis, and consist of "model- independent" and "model- 
dependent" terms. The model-independent terms are obtained from the nucleon-nucleon 
interaction, and by construction satisfy current conservation with it. The leading two- 
body weak vector current is the "7r-like" operator, obtained from the isospin-dependent spin- 
spin and tensor nucleon-nucleon interactions. The latter also generate an isovector "p- 
like" current, while additional isovector two-body currents arise from the isospin-independent 
and isospin-dependent central and momentum- dependent interactions. These currents are 
short-ranged, and numerically far less important than the 7r-like current. With the exception 
of the p-like current, they have been neglected in the present work. The model-dependent 
currents are purely transverse, and therefore cannot be directly linked to the underlying two- 
nucleon interaction. The present calculation includes the currents associated with excitation 
of A isobars which, however, are found to give a rather small contribution in weak-vector 
transitions, as compared to that due to the 7r-like current. The 7r-like and p-like contributions 
to the weak vector charge operator [15] have also been retained in the present study. 

The leading many-body terms in the axial current due to A-isobar excitation are treated 
non-perturbatively in the transition-correlation-operator (TCO) scheme, originally devel- 
oped in Ref. ||] and further extended in Ref. [26]. In the TCO scheme-essentially, a scaled- 
down approach to a full N+A coupled- channel treatment-the A degrees of freedom are 
explicitly included in the nuclear wave functions. The axial charge operator includes, in 
addition to A-excitation terms (which, however, are found to be unimportant [T5|]), the 
long-range pion-exchange term p7| , required by low-energy theorems and the partially- 
conserved-axial-current relation, as well as the (expected) leading short-range terms con- 
structed from the central and spin-orbit components of the nucleon-nucleon interaction, 



following a prescription due to Riska and collaborators [p8 



The largest model dependence is in the weak axial current. To minimize it, the poorly 
known iVA transition axial coupling constant g* A has been adjusted to reproduce the ex- 
perimental value of the Gamow- Teller matrix element in tritium /3-decay [§,0- While this 
procedure is model dependent, its actual model dependence is in fact very weak, as has been 
shown in Ref. [Bfl. 



The calculation proceeds in two steps [ 15[ : first, the matrix elements of p(q) and j(q) 
between the initial p 3 He LSJJ Z states and final 4 He are calculated with Monte Carlo integra- 
tion techniques; second, the contributing RMEs are extracted from these matrix elements, 
and the cross section is calculated by performing the integrations over the electron and 
neutrino momenta in Eq. ([I]) numerically, using Gauss points. 

The results for the S'-factor, defined as S(E) = E a(E) exp(4 tc a/v TC \) (a is the fine 
structure constant), at p 3 He cm. energies of 0, 5, and 10 keV are reported in Table |. In the 
table, the column labelled S includes both the 1 So and 3 Si channel contributions, although 
the former are at the level of a few parts in 10 3 . The energy dependence is rather weak: the 
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value at 10 keV is only about 4 % larger than that at keV. The P-wave capture states are 
found to be important, contributing about 40 % of the calculated S- factor. However, the 
contributions from D-wave channels are expected to be very small. We have verified explicitly 
that they are indeed small in 3 Dx capture. The many-body axial currents associated with 
A excitation play a crucial role in the (dominant) 3 Si capture, where they reduce the S- 
factor by more than a factor of four. Thus the destructive interference between the one- 



and many-body current contributions, first obtained in Ref. [Oj, is confirmed in the present 
study. The (suppressed) one-body contribution comes mostly from transitions involving the 
D-state components of the 3 He and 4 He wave functions, while the many-body contributions 
are predominantly due to transitions connecting the S-state in 3 He to the D-state in 4 He, or 
viceversa. 

It is important to stress the differences between the present and all previous studies. 
Apart from ignoring, or at least underestimating, the contribution due to P-waves, the 
latter only considered the long-wavelength form of the weak multipole operators, namely, 
their q=0 limit. In 3 P capture, for example, only the Co-multipole, associated with the 
weak axial charge, survives in this limit, and the corresponding S-i actor is calculated to be 
2.2 x 10~ 20 keV b, including two-body contributions. However, when the transition induced 
by the longitudinal component of the axial current (via the L -multipole, which vanishes 
at q=0) is also taken into account, the S-factor becomes 0.82 x 10~ 20 keV b, because of 
a destructive interference between the Cq and L RMEs. Thus use of the long- wavelength 
approximation in the calculation of the hep S'-factor leads to inaccurate results. 

Finally, besides the differences listed above, the present calculation also improves that of 
Ref. [|J in a number of other important respects: firstly, it uses accurate CHH wave functions, 
corresponding to the last generation of realistic interactions; secondly, the model for the 
nuclear weak current has been extended to include the axial charge as well as the vector 
charge and current operators. Thirdly, the one-body operators now take into account the 
1/m 2 relativistic corrections, which were previously neglected. In 3 Si capture, for example, 
these terms increase by 25 % the dominant (but suppressed) L\ and E\ RMEs calculated with 
the (lowest order) Gamow- Teller operator. These improvements in the treatment of the one- 



body axial current indirectly affect also the contributions of the A-excitation currents [|T5 
because of the procedure used to determine the coupling constant g* A . 

To conclude, we have carried out a realistic calculation of the hep reaction, predicting 
a value for the S'-factor five times larger than that used in the SSM. This enhancement, 
while very significant, is far smaller than that required by fits to the Super-Kamiokande 
data. Although the present result is inherently model dependent, it is unlikely the model 
dependence be so large to accomodate a drastic increase in the prediction obtained here. 

The authors wish to thank J.F. Beacom, J. Carlson, V.R. Pandharipande, D.O. Riska, 
P. Vogel, and R.B. Wiringa for useful discussions. The support of the U.S. Department of 
Energy under contract number DE-AC05-84ER40150 is gratefully acknowledged by L.E.M. 
and R.S. 
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TABLES 



TABLE I. The hep S-factor, in units of 10~ 20 keV b, calculated with CHH wave functions 
corresponding to the AV18/UIX Hamiltonian model, at p 3 He cm. energies £'=0, 5, and 10 keV. 
The rows labelled "one-body" and "full" list the contributions obtained by retaining the one-body 
only and both one- and many-body terms in the nuclear weak current. The contributions due 
to the S-wave channels only and S- and P-wave channels are listed separately. The Monte Carlo 
statistical error is at the 5% level on the total S- factor. 



E (keV) 











5 




10 






S 




S+P 


S 




S+P 


S 


S+P 


one-body 


26.4 




29.0 


25.9 




28.7 


26.2 


29.3 


full 


6.39 




9.64 


6.21 




9.70 


6.37 


10.1 
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